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^ I The status of the theoretical research on the compressional modes of finite nuclei and 

O ■ the incompressibility Koo of nuclear matter, is reviewed. It is argued that the recent 

■ experimental data on the Isoscalar Giant Monopole Resonance (ISGMR) allow extracting 
O ! the value of Koo with an uncertainity of about ± 12 MeV. Non-relativistic (Skyrme, 

Gogny) and relativistic mean field models predict for K^o values which are significantly 
different from one another, namely ~ 220-235 and ~ 250-270 MeV respectively. It is 
shown that the solution of this puzzle requires a better determination of the symmetry 

■ energy at, and around, saturation. The role played by the experimental data of the 
^ . Isoscalar Giant Dipole Resonance (ISGDR) is also discussed. 



1. INTRODUCTION 

The quest for the value of the nuclear incompressibility K^o is still continuing. Some 
significant progress in our understanding of how its value can be constrained have been 
achieved in recent times and this will constitute the subject of the present review. In this 
sense, this contribution is a continuation of those by J. P. Blaizot|l| and by N. Van Giai et 
a/.p] in the previous conferences of the Giant Resonance series (Varenna 1998 and Osaka 
2000). 

It is well known that the energy per particle E/A'm nuclear matter, considered as a 
function of the density exhibits a minimum at the saturation point Qq=Q.17 fm~^. The 
nuclear matter incompressibility, defined as 



(1) 



e=Qo 



provides a measure of the curvature oi E/A around Qq. The interest of determining K^o 
stems also from its impact on the physics of neutron stars. 

Since we cannot directly create and probe nuclear matter in ordinary laboratories, the 
only way to extract a value for Kao is by making contact with the phenomenology of 
the compressional modes in finite nuclei. The clearest example of compressional mode 
is the Isoscalar Giant Monopole Resonance (ISGMR), which is often called the nuclear 
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"breathing mode" and is excited by the operator 

M = j:rl (2) 

i=l 

The first evidences of this L = mode date back to the 1970s. The presence, in the 
same energy region, of modes with different multipolarities (e.g., L = 2) as well as of a 
non negligible background, makes the extraction of the monopole strength rather difficult. 
Over the last two decades, the experimental techniques have improved and our knowledge 
of the ISGMR properties has progressed. Many reactions have been employed to study 
this resonance, but inelastic [a, a') scattering has been, from the beginning, one of the best 
tools. In the light nuclei this resonance is rather fragmented, while in the medium-heavy 
nuclei it corresponds to a single peak of energy -Eisgmr ~ 80-A~^/^ MeV. Nowadays, in 
the recent measurements of ^°Zr, "^Sn, ^^^Sn and ^ospb performed at Texas A&M [Jj, the 
accuracy on the centroids of the ISGMR strength distribution has come down to about 
±2%. The importance of this high accuracy for the extraction of K^o will be discussed 
below. 



2. THE NUCLEAR INCOMPRESSIBILITY DEDUCED FROM THE IS- 
GMR: BASIC FORMULAS AND PREVIOUS RESULTS 

To perform the link between the properties of the ISGMR in finite nuclei and Koo, the 
definition (^Q) must be complemented by an operative expression which contains quantities 
that can be measured in the laboratory. J. P. Blaizot [A^ showed that a plausible definition 
of the finite nucleus incompressibility is given by 



- Vm<r2>' 

where m is the nucleon mass, and < > the ground state mean square radius. Using 
the experimental ISGMR energies to deduce for different values of A, there have 
been attempts to make an extrapolation of Ka to A=oo. This was done by using a 
Weiszacker-type formula for Ka, namely 

Ka = Koo + KsurfA-^^^ + KsyW + Kcoul^J^, (4) 

where a = M. Pearson [5^ was the first to show that, in view of the correlations 

among the parameters of (j3)) and the scarcity of experimental data, trying to use these 
to perform a fit of Eq. (@]) is statistically meaningless and would indeed leave Koo basi- 
cally undetermined (fits which lead to 100 MeV or 400 MeV may be equally acceptable). 
Similar conclusions were reached by S. Shlomo and D. Youngblood |6J. Therefore, we will 
not discuss these so-called "macroscopic approaches" to Koo. However, we remark that 
the theoretical values of the parameters entering Eq. can be calculated, within the 
framework of different models. We will come back to this point in Sec. 0] 

The procedure to obtain K^o, which is nowadays believed to be physically sound, is 
the so-called "microscopic approach" . The basic idea consists in using energy functionals 
E[g\ which allow calculating nuclear matter and finite nuclei on the same footing. 
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In the non-relativistic case, the starting point is a two-body effective nucleon-nucleon 
interaction V^g, whose parameters are adjusted to reproduce experimental data in a small 
set of nuclei. E is written as the expectation value of i^es = T + Ves on a independent 
particle wave function (i.e., a Slater determinant). In practice, the available functionals 
are based on the Skyrme and Gogny interactions. 

In the relativistic models, nucleons are described as Dirac particles which interact by 
the exchange of effective a, uj and p mesons. In the limit of large meson masses, point 
couphng models are obtained. In some cases, the coupling constants are taken to be 
density-dependent. As in the cases of Skyrme and Gogny, the couphng constants are 
fitted. A specific model has been recently proposed [7j, in which the main parameters are 
connected to more fundamental quantities, in particular to the so-called QCD sum rules 
and to the (iterated) pion exchange. In all models, the no-sea approximation is made and 
the energy functional is written in the Hartree (and not Hartree-Fock) form. 

In both non-relativistic and relativistic cases, the second derivative of the energy func- 
tional can be calculated analytically for uniform nuclear matter and the value of Koo 
associated to a given parametrization is therefore given. In the case of finite nuclei, one 
calculates the monopole excitation using self-consistent linear response theory. The system 
is perturbed with an (arbitrarily small) external field and the small oscillations around the 
ground state are governed by the residual force The theory is known as self-consistent 
Random Phase Approximation (RPA) and is well described in textbooks [H]. 

Then, the determination of K^o proceeds as follows. 

• Using a set of different parametrizations (within a given class of energy functionals) 
characterized by different values of ii'oo, self-consistent RPA calculations of the 
ISGMR are performed in a given nucleus. If the monopole strength has only one 
peak, -EisGMR is well defined and Eqs. (0111) suggest that a relation of the type 
-^iSGMR ~ \fKoo can be expected. This can be verified empirically and indeed, 
relations of the type 



have been interpolated (see below). 
• The experimental value of -Eisgmr is inserted in Eq. ^ and the value of K^o is 



One or few nuclei would be enough to apply this procedure, and ^ospb jg a typical system 
where the monopole strength has a well-defined peak and where recent experiments have 
reduced the sources of errors. From the theoretical point of view, the calculations must 
be free of the uncertainites associated to, e.g., the description of pairing or to anharmonic 
effects. 

The procedure that we have described, was firstly applied by J. P. Blaizot and collab- 
orators jH], by employing the Gogny interaction. The authors of made use of the 
existing parametrizations and also built ad hoc new ones in order to cover more values 
of Koo- A value of about 230 MeV for the nuclear incompressibility can be extracted 
from the experimental -Eisgmr of ^°®Pb. Exactly the same procedure was applied in the 




(5) 



deduced. 
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case of the Skyrme forces (using only already existing parametrizations) jim2j . From the 
experimental monopole energy of ^°*Pb a value of K^o around 210 MeV is deduced, while 
using ^°Zr the value is even lower (around 200 MeV). All the RPA calculations quoted 
are done using a discrete basis. However, the value of 210 MeV is consistent with what 
has been found by other authors within the framework of continuum- RPA [TT] . 

Different calculations made by using the relativistic RPA gave instead larger values of 
Koo- Values of 250 MeV and 270 MeV were extracted from the experimental data in 208p|-j 
and ^^"^Sm respectively |12|13j . This model dependence in the extraction of K^o has been, 
and partly it is still, the basic puzzle. 

3. CONSISTENCY BETWEEN THE RESULTS FROM SKYRME AND GOGNY 
INTERACTIONS 

Before proceeding, some quantitative statements concerning the numerical accuracy of 
the microscopic approach are in order. From Eq. (0), the relative errors on the ISGMR 
energy and on Koo are related by 

5Kao _ 2^^ISMGR 
Koo -f'ISGMR 

As a rule of thumb, let us keep in mind that sticking to the case of ^''^Pb, ±150 keV of 
uncertainity on -Eismgr result in about ±5 MeV uncertainity on K^o- The experimental 
measurement on 208p]^ provides us with 14.17 ± 0.28 MeV fS^ and therefore 

5KZ^- ~ ±10 MeV. 

Theoretically, the best way to extract the monopole energy is by means of constrained 
Hartree-Fock (CHF) calculations ^1]. These calculations provide the inverse energy- 
weighted sum rule m_i with a numerical error of the order of ±3% (see also Ref. |15j). 
and since the energy-weighted sum rule mi is known, the relation 



-ElSGMR — J — ^ (6) 

implies 

~ ±7 MeV. 

The two errors on K^o are independent and should be added quadratically, so that 

5Koo ~ ±12 MeV. (7) 

We now discuss in more detail the calculations performed in Refs. jl0|2j . One of their 
main limitations is the lack of complete self-consistency, as the residual p-h Coulomb and 
p-h spin-orbit interactions are dropped. We have analyzed the effect of this approximation 
in the nuclei ^^O and ^°Ca. The results are shown in Tabled The Skyrme force used is 
SLy4 jni- The results of (a) are obtained by dropping the Coulomb and spin-orbit terms 
both in HF and in RPA. In this sense, the calculation is fully self-consistent, and the fact 
that the RPA results agree with the CHF results within less than 1% suggests that the 
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Table 1 

Values of the m_i sum rule in fm^/MeV calculated using the interaction SLy4. In column 
(a), results obtained without the Coulomb and spin-orbit interaction are shown. In this 
case, the RPA is fully self-consistent. In column (b), the Coulomb and spin-orbit terms of 
the interaction are included in the CHF calculation and in the HF calculation on which 
RPA is based, but are excluded from the residual RPA interaction. 







(a) 


(b) 


160 


CHF 


14.45 


16.04 




RPA 


14.49 


16.73 


^oCa 


CHF 


75.31 


88.31 




RPA 


75.91 


92.13 



approximations done in the RPA (neglect of the continuum and truncation of the discrete 
basis) do not affect seriously the values of m_i. On the other hand, the results of (b) are 
obtained by including the Coulomb and spin-orbit terms in the mean field and not in the 
residual RPA interaction. The difference between the CHF and RPA results is larger than 
the CHF intrinsic uncertainity and can therefore be a meaningful indication of the error 
induced by the lack of full RPA self-consistency. Consequently, we have compared CHF 
and RPA in the nuclei already considered for the extraction of K^oi that is, ^°^Pb and 
^°Zr. Note that in the case of ^osp^ ^^Xt the CHF result for m_i is larger than the 
RPA result, contrarily to the case of ^^O and ^°Ca. The results are displayed in Figs. ^ 
and 121 

The values of 235 MeV from ^osp]^ g^j^^ 220 MeV from ^°Zr for K^o which are extracted 
from the Skyrme-CHF calculations, are consistent with each other and with the values 
deduced using the Gogny interaction. This is the first important conclusion of the present 
paper. 

We conclude this Section with a brief justification of the mean field approximation in 
the present calculations. Firstly, it has been shown in Ref. that including the coupling 
of the RPA states with more complicated configurations of 2 particle-2 hole (2p-2h) type, 
shifts the ISGMR in ^'^^Pb to lower energy by only ^ 500 keV. This number is smaller 
than the uncertainities discussed above. Secondly, according to the arguments of Sec. |2l 
in order to use the result for Ejsgmr obtained beyond mean field, a calculation for nuclear 
matter performed on the same footing is needed. The development of a suitable energy 
functional is still a challenge for nuclear structure theory. 

4. HOW TO RECONCILE THE RELATIVISTIC MEAN FIELD WITH SKYRME 
AND GOGNY ? 

As far as the relativistic models are concerned, there has been a recent suggestion jTH| 
that the different outcome for (compared with Skyrme and Gogny) is originated by the 
different density dependence of the symmetry energy S{p) predicted by different models. 
It is true that the measurements of the ISGMR are done in a system (^'^'^Pb) with a finite 
value of (N — Z)/A. To illustrate the possible influence of the density dependence of 
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Figure 1. Monopole energies in ^°^Pb (defined as in Eq. ©) obtained from RPA (crosses) 
and CHF (stars) calculations wliicli employ different Skyrme forces, plotted as a function 
of the associated K^. The lines are fits of the type (0). The dashed lines indicate the 
extracted values of K^o- The value of K^o resulting from CHF is in agreement with that 
extracted from the Gogny calculations of 

S{p) on the extraction of i^oo, J- Piekarewicz ^Hj has built parametrizations of effective 
Lagrangians whose symmetry energy has different density dependences (this is easy to 
achieve since the p coupling constant is an adjustable parameter) and he finds that the 
extracted K^o indeed differ and can even become close to Skyrme force values. To achieve 
this, it is necessary to soften the function S{p) and in Ref. 18] this was done by lowering 
the symmetry energy at the saturation point, a,-. However, it has been pointed out in 
Ref. that, in the type of model used in Ref. parametrizations with a,- lower than 
36 MeV cannot describe satisfactorily N ^ Z nuclei. A complementary attempt has been 
made in [15], by constructing Skyrme forces with associated values of ar up to 38 MeV. 
In this case, finite nuclei have been carefully considered in the fits. By calculating the 
ISGMR centroid energy, a very weak dependence on ar has been found. 

There is one more conceptual remark. Starting from (j3)), one would expect a dependence 
of Kp^, and consequently of the monopole energy, on a parameter like Kgym more than 
on a-r- That is, on the derivatives of the symmetry energy more than on its value at 
saturation. In fact, it has been found that both in relativistic and in non-relativistic 
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Figure 2. The same as Fig. ^in the case of ^°Zr. 



models the quantity Kg^rf is essentially given by cKoo with c ~ -1 |4|2Uj . Therefore, we 
can write 



^ sym ^ ' Coul A^/"^ ' 

i^A ~ i^rni + cA-V3) + i^(-^)«2 + 4-;;)^. (8) 

It is likely that the third term of the r.h.s. (Coulomb contribution) does not change much 
from a non-relativistic to a relativistic description. The same values of can thus be 
obtained with different values of the and Kgym terms. For illustration, we display in 
Fig. El the results of the Skyrme calculations of the ISGMR in ^°^Pb (already shown in 
Fig. 1) together with the corresponding results of the relativistic mean field taken from 
|12j . In this nucleus we have = 0.04. At any given value of Koo the difference between 
the two curves is approximately 1 MeV. This translates into a difference of about 20 MeV 
in Ka- If this difference is entirely attributed to the negative term KgymOp' of Eq. (jSJ, 
the values of Ksym in the non-relativistic and relativistic models would differ by about 
500 MeV. Only few calculations of Ksym are available |4f2()j . therefore more systematic 
tests of the present argument should be made. More importantly, we apparently miss any 
experimental constraint to decide what is the proper value of the parameter Kgym- 
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Figure 3. Monopole energies as a function of in Skyrme and in relativistic mean 
field (RMF) models. The relativistic results are taken from Ref. jTSj. The Skyrme results 
are the same as in Fig. ^ (only the point corresponding to the SIII interaction has been 
added in order to compare with the relativistic result for K^o ~ 350 MeV). The lines are 
numerical fits. 



5. THE ISOSCALAR GIANT DIPOLE RESONANCE 

The ISGDR is a non-isotropic compressional mode, which is excited by the operator 

b = Y.rff^M{h). (9) 

i=\ 

and which provides in principle a further way to extract the value of . There are many 
reasons why this task is more challenging in the case of the dipole than in the case of the 
monopole. 

From the theoretical point of view, one difficulty arises from the fact that the spurious 
center-of-mass translation (associated with the operator Yld=\'f'i^\Miji)) carries the same 
quantum numbers as the ISGDR. In principle, a sharp spurious state at zero energy should 
result from an ideal self-consistent RPA calculation. However, in practice this is not the 
case due to approximations and numerical inaccuracy. Consequently, the resulting states 
are not orthogonal to the spurious state and one has to correct for this. The spurious 
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transition density is expected to be of the type ~ -Jr where is the ground state density. 
It is possible to project out the spurious component from each excited state (cf., e.g., 
Ref. [21] )• An equivalent procedure [22] consists in using instead the modified operator 

A 

Deff = T.irf-Vri)YiMiri), 

i=l 

where = | < >. A recent analysis of the accuracy of these techniques, as well 
as a discussion about some different prescriptions to subtract the spurious state, can be 
found in Ref. [2S1- There have been other discussions in the recent literature [21], also in 
connection with semiclassical models 

On the experimental side, different {a, a') measurements have been performed over the 
years but the problem of disentangling the ISGDR strength from the other multipoles (and 
from the IVGDR) is far from being trivial. The ISGDR lies at higher energies than the 
ISGMR (approximately 110-A~^/^ MeV). An accurate determination of its high-energy 
tail is therefore more difficult. On the low-energy side, a sizeable amount of fragmented 
strength is found. These issues are discussed in these Conference proceedings [2612"?] . 

Different theoretical calculations |21l28j . have clarified that the low-energy part of the 
ISGDR strength is formed by non-collective states. The separation of a "high-energy" 
region and a "low-energy" region in the ISGDR strength distribution, emerges system- 
atically from the calculations of Ref. [2I| as illustrated here in Fig. HJ One of the main 
indications about the different character of the two parts of the strength, comes from the 
fact that the centroids of the high-energy regions, calculated with different Skyrme forces 
in a given nucleus, scale with the corresponding Koo (which testifies to their compres- 
sional nature), whereas the centroids of the low-energy regions do not. The same pattern 
is found in the relativistic calculations. 

More detailed considerations concerning the wave functions of the low-lying states have 
been done in Ref. [SH] where the authors suggest that the low-lying strength may corre- 
spond to a "toroidal" resonance, which can be visualized - in a simple way - by thinking 
of a particle current bent into a torus. It is excited by the vector operator 

A 

f = ^ V X (rl X V)rf FiAf(f,) (10) 
1=1 

which couples to the transition current. To verify the toroidal nature of the states, some 
other probe than the a-particles should be tried, and (e, e') experiments would be probably 
useful in this respect. 

In Table 121 we report the results of the various self-consistent RPA calculations of the 
ISGDR in ^°^Pb. The first remark is that the low-energy part is of course significantly 
dependent on the model and on the specific functional used, as expected due to its lack 
of collectivity. As far as the high-lying centroid is concerned, the Skyrme results of 
Refs. |3H21j (continuum and discrete RPA respectively) are in good agreement with each 
other. In Ref. [HI] it has been shown that coupling the RPA states with 2p-2h type 
configurations, shifts the ISGDR centroid down to 22.9 MeV, in very good agreement 
with experiment (this coupling also produces a conspicous spreading width of about 6 
MeV). 
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Figure 4. ISGDR strength functions in different nuclei, calculated in RPA using the 
Skyrme interaction SGII 29] and corrected for center-of-mass effects. In the case of ^osp]-, 
the dashed line corresponds to a calculation where the spurious center-of-mass state is 
not subtracted. Here as well as in Fig. 5 the discrete RPA states have been smeared out 
by means of a Lorentzian of 1 MeV width. Taken from Ref. ||21j. 



In Table 2, also the results of the most recent measurement performed at RCNP in 
Osaka (making use of a-particles at incident energy of 400 MeV) as well as the new 
findings of the Texas A& M group, are shown. The experimental energies seem to converge 
towards each other (compared to the previous experiments of the same groups, where 
the discrepancies were larger). Nevertheless, a point should be made concerning the 
experimental analysis. The directly measured quantity is the double-differential cross 
section, J^^^ - The multipole decomposition of this cross section is done by relying on 
Distorted Wave Born Approximation (DWBA) calculations where the radial form factors 
for the various multipole excitations are the same at all energies. This is somehow in 
contrast with the outcome of theory, as we have discussed above. 

We conclude the discussion on ^ospb with a statement about the deduction of Koo- If 
this quantity is derived from a plot similar to those of Figs. [T] and El i.e., by using the 
"microscopic approach" described in Sec. Elin connection with the ISGMR but using now 
the ISGDR data, one finds a value around 205 MeV which is still (marginally) compatible 
with the value 220-235 MeV discussed above (cf. Fig. 2 of [21j). This means that there 
is no basic contradiction between the present ISGDR and ISGMR data, as far as the 
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Table 2 

Self-consistent (relativistic and non-relativistic) RPA calculations performed for the IS- 
GDR in ^''^Pb compared with the most recent experimental data. The two columns report 
the centroid energies (in MeV) of the low-energy and high-energy regions discussed in the 
text. 

High-energy Low-energy 

centroid centroid 
Hamamoto et al. |,31j 23.4 ~ 14 

Colo et al. |21j 23.9 10.9 

Vretenar et al. [2HI 26 10.4 

Piekarewicz |S2! 24.4 ~ 8 

Shlomo and Sanzhur EH ~ 25 ~ 15 



Uchida et al. [33-2H] 23 ± 0.3 12.7 ± 0.2 



Lui et al. EH 21.7 12.6 



deduction of the nuclear incompressibility is concerned. However, it should be clear from 
all our discussion, that we still need to await for improvements in the ISGDR studies in 
order to reach the same confidence that we have in Koo extracted from the ISGMR. 

We end this Section by showing the result of a Skyrme-RPA calculation for a lighter 
system, namely ^^Ni (see Fig. E)). It can be noticed, by comparison with Fig. 01 that 
although the centroid of the ISGDR is still in the region around llO-A-^/^ MeV, the 
strength is more fragmented than in the heavier systems. Also, the distinction between 
a low- and a high-energy regions is less evident. A comparison with experiment would 
be useful. In any case, as it was said for the monopole case, this picture confirms that 
medium-mass systems are less suitable for studying the nuclear incompressibility. 

6. CONCLUSIONS 

In recent years, there have been significant progresses both in the experimental tech- 
niques aimed to extract with good precision the moments of the monopole strength func- 
tion, and in the theoretical models - especially in those based on relativistic functionals 
which can be nowadays discussed on the same footing as the well tested Skyrme and Gogny 
functionals. Still, it is shown in this paper that it is in general possible to determine a 
delicate quantity like the nuclear incompressibility only within ±12 MeV. 

Many discussions have been devoted to the fact that the values of Koo depend on the 
model through which they have been extracted, either Skyrme, Gogny or relativistic func- 
tionals. In this contribution, we show that the discrepancy between Skyrme and Gogny 
does not exist and, using the ISGMR data in ^°^Pb the value of K^o lies between 220 MeV 
and 235 MeV. On the other hand, the relativistic calculations point to larger values, of 
the order of 250-270 MeV. This puzzle is still unsolved, but we argue in this work that the 
reason has to be found probably in the different features that the asymmetry energy curve 
has in the relativistic and non-relativistic models, respectively. The relativistic models 
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Figure 5. ISGDR Strength distribution in ^^Ni calculated using Skyrme-RPA with the 
SGII force. 



are characterized by significantly larger values of the symmetry energy at saturation, and 
of its first and second derivatives. The quest about the proper values of the asymmetry 
energy has of course very general implications on the nuclear phenomenology, and new 
works on the subject are in progress. 

As far as the isoscalar dipole is concerned, there has been also definite experimental 
improvements in the techniques aimed to reduce or eliminate the background, and in the 
data analysis. While in the past the results from different groups were in disagreement 
among themselves, and with theory, this does not seem to be the case for the new ex- 
periments. The degree of accuracy that can be expected in extracting K^o is smaller in 
the dipole than in the monopole case. However, there does not seem to exist a basic 
incompatibility between the values of the incompressibility deduced using for instance 
the Skyrme forces, in the monopole and in the dipole case. Among the perspectives, it 
can be said that new {a, a') cross section calculations based on a fully microscopic input 
may be of great help to reduce some of the uncertainities which still exist. On the other 
hand, it is unlikely that the ISGDR data can help in solving the discrepancy between the 
non-relativistic and relativistic models as far as the value of Koo is concerned. 
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